Dynamic freezing of strongly correlated ultracold bosons 
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We study the non-equilibrium dynamics of ultracold bosons in an optical lattice with a time 
dependent hopping amplitude J{t) — Jq + SJ cos{ujt) which takes the system from a superfluid phase 
near the Mott-superfluid transition (J — Jo + 5 J) to a Mott phase {J — Jo — 5 J) and back through 
a quantum critical point (J = Jc) and demonstrate dynamic freezing of the boson wavefunction 
at specific values of oj. At these values, the wavefunction overlap F (defect density P = 1 — F) 
approaches unity (zero). We provide a qualitative explanation of the freezing phenomenon, show 
it's robustness against quantum fluctuations and the presence of a trap, compute residual energy 
and superfluid order parameter for such dynamics, and suggest experiments to test our theory. 
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Theoretical study of non-equilibrium dynamics in 
closed quantum systems has seen great progress in re- 
cent years ^ mainly due to the possibility of realiza- 
tion of such dynamics using ultracold atom in optical 
lattices [HIS]- For bosonic atoms, such systems are well 
described by the Bose-Hubbard model with on-site inter- 
action strength U and nearest neighbor hopping ampli- 
tude J m [S] . Several theoretical studies have been car- 
ried out on the quench and ramp dynamics of this model 
[6llllj ; some of them have also received support from re- 
cent experiments 0. In contrast, studies on periodically 
driven closed quantum systems have been undertaken in 
the past mainly on driven two-level systems [I2l [13] or 
on weakly interacting or integrable many-body systems 
which can be modeled by them [Ml [15]. Among these, 
Ref. T5] has predicted freezing of the time-averaged value 
of the order parameter (magnetization) of an periodically 
driven one-dimensional (ID) Ising or XY model, when 
the temporal average is performed over several drive cy- 
cles, at specific drive frequencies. Such a freezing occurs 
in the high frequency regime and exhibits non-monotonic 
dependence on the drive frequency. However, to the best 
of our knowledge, the phenomenon of dynamic freezing 
has never been demonstrated for dynamics involving a 
single drive cycle and/or for non- integrable quantum sys- 
tems. Recent studies of periodic dynamics of the Bose- 
Hubbard model have not addressed this issue [TBI E] • 

In this work, we demonstrate, via designing a peri- 
odic driving protocol, that the periodically driven Bose- 
Hubbard model may exhibit dynamic freezing of the bo- 
son wavefunction \tfj{t = 0)) = \tp(t = T)) for specific 
values of the drive frequencies lu = 2tt/T. Our driv- 
ing protocol constitutes a time-dependent hopping am- 
plitude of the bosons J{t) — Jq + SJcos{ujt) with Jq 
and 6 J chosen such that the drive takes the system from 
a superfluid (SF) (J — Jq + SJ) to the Mott insulator 
(MI) state (J = Jo — SJ) and back through the tip of 
the Mott lobe where fi — ^tip- We demonstrate, using 
mean-field theory, that such a freezing phenomenon de- 



rives from quantum interference of the dynamic phases 
acquired by the bosons and compute the defect formation 
probability P = 1-F (where F = |(-0(i = 0)\^p{t = T))\'^ 
is the wavefunction overlap), the superfluid order pa- 
rameter A(r) = {ip{T)\b\^l;{T)) (where b denotes the 
boson annihiliation operator), and the residual energy 
Q^E{t = T)- Eg (where E{t = T) is the energy of the 
system at the end of the drive cycle and Eq is the initial 
ground state energy) as a function of uj. We also show, 
via inclusion of quantum fluctuation by a projection op- 
erator approach jllj and numerical mean-fleld study of 
a trapped boson system that the freezing phenomenon 
is qualitatively robust against quantum fluctuations and 
the presence of a trap. We note that such a freezing be- 
havior has two novel characteristics which distinguishes 
it from its counterpart in Ref. [T^ . First it does not need 
high frequencies a,s huj/U ^1 throughout the range of 
uj where the freezing occurs. Second, it occurs for a sin- 
gle cycle of the drive and does not need averaging over 
several cycles. Such a dynamic freezing phenomenon has 
not been studied in the context of closed quantum sys- 
tems; our work therefore constitutes a significant advance 
in our understanding of periodic dynamics of closed non- 
integrable quantum systems. 

The Hamiltonian describing a system of ultracold 
bosonic atoms confined by a trap and in an optical lattice 
is given by 

n = -Jb%, + ^[-Airrir + ^nr(nr - 1)] (1) 

(r,r') r 

where denotes the chemical potential at site r, r' de- 
notes one of the z nearest neighboring sites of r, and 
fij. — b\br- In the absence of a trap, fi^ = n for all sites 
and for zJ <^ U, the ground state of the model is a MI 
state with n bosons per site with n — 1 for < fi/U < 1. 
For zJ 3> U, the bosons are delocalized and the system, 
for d > 2, is in a SF state. In between, at J = Jc, the 
system undergoes a SF-MI transition. The equilibrium 
phase diagram of the model constitutes the well-known 
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FIG. 1: (Color online) Top left panel: Plot of 1 - ri(t) (red 
dashed line) and (/>s(f) (blue solid line) as a function of t. 
Bottom left panel: Variation of 1— ri(T) (red dashed line) and 
4>b{T) (blue solid line) with lj. Top (Bottom) right panel: Plot 
4>dit) {4>d{T)) as a function of t (cj). For all plots Jo — 1.05Jc, 
5J = 0.35Jc, and = 0.414C/. 



Mott lobe structure [H [S] . 

To obtain an semi-analytic insight to the freezing phe- 
nomenon, we first analyze the periodically driven Bose- 
Hubbard model in the absence of a trap and within 
mean-field approximation. The time-dependent mean- 
field Hamiltonian is given by 

Hnif = ^[-Ainr + ^nr(nr-l)] + (A;(<)&t +h.cX?) 



where A;(t) = -J{t)J2{r'){f^r') and Ao = A'(t = 0). 
Within homogeneous mean-field theory, the Gutzwiller 
wavefunction for the bosons reads |'0(r,t)),„f — 
nrSn/"(^)l^) [IS|- The Schrodinger equation 
idt\ip{r,t))mf = ?^mf(OIV'(r,i))mf yields the time- 
dependent mean-field equations for /„(i) = /„: 



{idt - E,,)f„ = A(^)^^/„_l + A*it)Vn + 1/„h 



(3) 



where A{t) = -zJ{t)J2nV^fn-ifn, and En = -nn + 
U{n~ 1)71/2. In what follows, we shall choose Jo and 5J 
such that the ground state of J/mf with J = Jo -I- ^ J is 
a SF state close to the QCP so that /„(t = 0) ~ for 
n > 3 and fi{t = 0) > fo{t = 0) ,/2(i = 0). As can 
be verified by explicit numerics, in this regime /n(t) for 
n > 3 remains small for all t during the dynamics and can 
thus be neglected. The equations for /o, /i and /2 then 
reads (suppressing time dependence of /„(i) for clarity) 



idth 
idtf2 
idth 



= -zJit)[\h\'fo + V2f;fi] 

= E,h-zJit)[2\h\y, + V2f*f^] 



(4) 




FIG. 2: (Color online) Top left panel: Plot of the defect den- 
sity P as a function of u/Jc displaying freezing at a; = cj^. 
Inset shows Logj^gP near ujg ~ 0.47Jc. Bottom left panel: 
Variation of Urn/Jc as a function of Jo/Jc for m — 6. Right 
panels: Plot of Q and |A(r)| as a function of u/Jc. All pa- 
rameters are same as Fig. [T] 



From Eqs. [4j it is easy to see that |/„|, for n < 2, obeys 



the relation dt\fo\ - dt\f2\ 



-9t 1/1 1/2. Parameteriz- 



ing /„ — rn{t) exp[i(/)„(t)], with the choice that (/)„(0) = 
and ^^(O) = fn{0), one can write 



'2[0] 



it) 



-(r?(t)-l)/2+H^, 



(5) 



where 77 is a time independent parameter whose value is 
fixed by the initial values r„ [12]. Note that rj represents 
the magnitude of the particle-hole asymmetry since Tq = 
r2 for 77 = 0. Substituting Eq. [5] in Eq. |4j we get 

dtri = -V2zJ{t)sin{(l)s)rigo{ri), 

dtcjjs = -U + zJ{t)[gi{ri) ~ g2{ri)cos{(j3s)\-, (6) 



U + 2^l + zJ{t)rl 1 - 4y27ycos((/)^)/go(fi) 



where we have suppressed the time dependence of ri and 
for clarity, = (?!)o + '/'2 - 2</)i and (jy^ = (j)2 - 
(po are the sum and differences of the relative phases of 
the Gutzwiller wavefunction, and the functions gi(ri) are 
given by 



go(ri) - y^(l-r2)2 - 4772, g,{n) = Grj - 3 - 2fj, 

32 (ri) - 2V2 [rlirl - l)/<?o(ri) + 50(^1)] . (7) 

We note that the first two of the equations in Eq. |6] are 
coupled equations describing the evolution of ri and 0s, 
while the third describes the evolution of (/)ci in terms 
of ri and (jfs. Furthermore, using a scaled variable f — 
wt/(27r), we find that the relation between ri and (j)s can 
be written as 



E,h - zJ{t)\{2\h? + l/oH/i + 2\/2/*/2/o]. 



dri/d(t)s 



^\/2sin{(j)s)rigo{ri) 



[gi(ri) - .g2(ri)cos(0,)] - U/zJ{t') 



TV (8) 
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FIG. 3: (Color online) Similar plots as in Fig. [5] but with 
Jo[SJ] — 0. 87[0. 58] Jc and computed using the projection op- 
erator approach displaying robustness of the freezing phe- 
nomenon against quantum fluctuations. See text for details. 



This allows us to symbolically write — ^{Ti,t'), where 
^ is an unknown function, and thus establish an uj inde- 
pendent relation between r and 0s for any fixed t' . 

Eqs. [5j [6] and [7] constitute the central result of this 
work. They constitute a complete description of the evo- 
lution of /o, /i, and /2 in the presence of the periodic 
drive and provide an understanding of the freezing phe- 
nomenon as follows. First, we find that a numerical solu- 
tion of Eq. [6] together with Eq. [5] allows us to obtain r„, 
and as a function of time. A plot of 1 — ri{t) and 
4>s{t) as a function of t for lu/Jc = 0.52 is shown in the left 
panel of Fig. [l] We find that ri changes appreciably when 
J{t) is close to Jc, however, ri(T) ~ fi(0) at the end of 
the evolution. Note that this also implies, via Eqj5j that 
^2(0) (0) — ^2(0) (^)- The bottom left panel of Fig.lllshows 
that the relation ri{T) ~ ri{0) holds for a significant 
range uj/Jc < 0.8. Second, we note that, 4>s{t) undergoes 
rapid oscillation when J{t) < Jc] however, it also comes 
back close to it's initial value at the end of the drive: 
4>s(T) ~ f^sCO). Since (j)s and ri satisfies a oj independent 
relation, (j)s = ^(ri,i'), we infer that 0^ must remain 
close to its initial value for the same range of uj for which 
ri(T) ~ ?'i(0); this is verified numerically in the bottom 
left panel of Fig. [T] Finally, we note, from the right panels 
of Fig. [l] that 4>d{T) is a monotonic function of uj. Thus 
we may define w = for which 4>d{T) — 47rm (m being 
an integer). Together with the fact that ri(T) ~ ?'i(0) 
and (t>siT) ~ (/>s(0) — 0, we find that at uj = both the 
relative phases satisfy (t>2 — (f'l = ~i4'o ~ — Zttto lead- 
ing to \^Jnif{T)) ~ |^/'mf(0)) up to a global phase. This 
constitutes the dynamics freezing of \4')mf- 

To obtain an accurate estimate of the degree of freez- 
ing, we compute the defect density P = 1 — F = 
1 - l(0mf(7')IV'mf(O))p. The plot of P as a function of uj 
clearly shows that P — )• at w = w^. A plot of Log^gP 
vs UJ near ujg/ Jc ~ 0.47, shown in the top left panel of 



Fig. [2] reveals that P ^ 10~^ indicating that the over- 
lap, up to a global phase, is exact within our numerical 
accuracy. We have checked for all uj^ < 0.8 Jc, P < 10^'^ 
which indicates a near perfect freezing. We also compute 
the residual energy Q{T) and the SF order parameter 



A 



(9) 



at t = T as a function of uj. We find from Eq. |9] that 
|A| is independent of 0^. Thus |A(T)|/|Ao| and Q/U 
(which can also be shown to be independent of (jjd) remain 
close to unity and zero respectively over the entire range 
of uj/ J for which ri and (j)^ remain close to their initial 
values as shown in right panels of Fig. [2j Such a behavior 
distinguishes these quantities from P which depends on 
(j)d and hence vanishes at discrete uj„i. Finally, we find 
that for all values of Jq shown in bottom left panel of 
Fig. [2] there is an appreciable range ofuj/Jc within which 
the freezing phenomenon occurs and that decreases 
monotonically as a function of Jq over this range. 

Next, we study the effect of quantum fiuctuations on 
the freezing phenomenon. We incorporate such fluctua- 
tions by using a projection operator method developed 
in Ref. JJj which provides an accurate treatment of dy- 
namics with fluctuations for zJ{t)/U ^ 1. The idea 
behind this approach, as detailed in Ref. [TTJ is to intro- 
duce a projection operator Pi — |?T.o)(no|r x I'^o) ("o|r' 
which lives on the link i between the neighboring sites 
r and r' of the lattice. Using Pg, one can write the bo- 
son hopping term as T' = X](rr') -J{t)blbr' = J^i^i = 

Eemn + nPi) + PtnPh where Pt = (1 - Pt). 
In the strong-coupling regime where zJ{t)/U ^ 1, the 
term Tg^[J] — [PiT'g + T^Pe), at any instant, represents 
hopping processes which takes the system out of the 
instantaneous low-energy manifold. Thus one can de- 
vise a time-dependent canonical transformation via an 
operator S = S[J{t)] = Y.i-APi^T'A/U which elimi- 
nates Ti^[J{t)] up to first order in J{t)/U and leads to 
the effective instantaneous time-dependent Hamiltonian 
H* ~ exp{—iS[J{t)])Hexp{iS[J{t)]). Such a canonical 
transformation is equivalent to a transformation on the 
system wavefunction \ip): \ijj') — exp{-'iS[J{t)])\tp) . We 
note that and 1-0') coincides for J = which leads 
us to the natural choice lip') = Hr X^n /r^lOI^-r)- Note 
that \ip) is not of Gutzwiller form; it involves spatial cor- 
relation due to exp{iS[J{t)]) factor. The instantaneous 
energy of the system is given by _E[{/"(i)}] — {^p\'H\'^^}) — 
{iP'\H*\ij') + 0{J{tf/U^) and includes 0{J^/U^) quan- 
tum fluctuation corrections. As shown in Ref. 11 , this 
formalism allows one to describe the dynamics of the 
bosons by solving for the Schrodinger equation for \ip'): 



{ihdt + dS[J{t)]/dt)\i;') = H*[J{t)]\i>'} 



(10) 



Using the expression of \iJj') and i?[{/„}], one can convert 
Eq. 10 to a set of equations for {/"(t)} pi]. Defining 
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FIG. 4: (Color online) Top left panel: Plot of the boson den- 
sity profile in the trap at t = and t = T. Top right panel: 
Variation of relative phase cf>d{T) as a function of the position 
(along 2/ = 0) in the trap displaying coherent evolution of the 
bosons. Bottom panels: Plot P and |A(r)|/|Ao| as a func- 
tion ui displaying freezing phenomenon at ui^. For all plots 
Jo = 0.04i7, 5J = 0.015(7, Jc = 0.041(7, and /io = 0.415(7. 



fpn = \Jn + l/*/„+i, one gets 

ihddu = <5i?[{/„(t)};J(t)]/(5/: 



izh dJ{t) 

1j dT 



A numerical solution of Eq. 11 yields fn{t) and hence 
IV'') using which one can compute \'4'[t)) = expfiS*] |V''(t)) 
perturbatively to 0[J(t)^/J7^]. Similarly, expectation 
value of any operator O at any instant t can be calcu- 
lated in terms of W{t)): (O) = {ij'{t)\e-'^Oe'^\ip'it)) = 
{■ilj'it)\0\^p'{t))-{ilj'{t)\[iS,0]\ip'it)) + ..., where the ellip- 
sis indicate higher order terms in J{t)/U. Note that the 
second term in the expression originates from quantum 
fluctuation and modifies mean-field result (first term). 
Using the above-mentioned procedure detailed in Ref. 
dH, we compute P(T), Q{T) and |A(T)| as shown in 
Fig. [3) We find that key effects of the quantum fluctua- 
tions is to change numerical values of wjjj and the precise 
range of w over which freezing occurs; however the mean- 
field results hold qualitatively in the sense that P — ^ 
for several with Log^pP < —4 for all w^. Further 
o;^ also decreases monotonically with Jq as shown in left 
bottom panel of Fig. [3] for uj* ~ 0.6Jc. 

Finally, we consider the effect of a harmonic trap on 
the freezing phenomenon. For this part, we numerically 
solve Eq.[3]for d^2 with /^r = i-i-o + Q.OlU[{v^ - l/2f + 
(ry — 1/2) ], for Nq = 576 sites (linear dimension 24) and 
with fixed total particle number Nq. We choose the trap 
parameters so that the ground state of the bosons in cen- 
ter of the trap at i = r/2 is MI phase with n = 1. The 
evolution of the density profile of the bosons is shown in 



the top left panel of Fig. [I] for t = (left) and t = T/2 
(right). The top right panel indicates evolution of 0^ as 
a function of the position of the bosons in the trap along 
the line y = 0. The plot indicates that for all uj < Jc, 
(j)d evolves coherently with negligible spatial variation. 
The plots for (ps and ri are similar in nature; thus, we 
expect the boson evolution to have the same qualitative 
properties as that found within a homogeneous mean- 
field approach. A plot of P (|A(r)|/|Ao|) as a function 
of Log2o(w/ Jc) in the lower left (right) panels of Fig. [4] 
confirms this expectation. We find that the main effect of 
the trap is to push the freezing phenomenon to lower fre- 
quencies leaving its qualitative nature unchanged. The 
largest freezing frequency occurs at ~ 0.2 Jc which is large 
compared to frequencies ~ 0.05 Jc where momentum con- 
serving boson pair production at finite momenta, which 
is not captured within mean-field theory, is expected to 
become significant [2^. Fig. |4] also demonstrates that the 
freezing phenomenon disappears at higher drive frequen- 
cies where the trapped bosons do not evolve coherently 
leading to spatial variation of 4>d- 

For experimental verifications of our work, we suggest 
interference of two bosonic condensates in the presence of 
an optical lattice, near the QCP which are separated after 
creation by a double-well potential and allowed to evolve 
separately for a fixed holdout time. It is well known that 
recombination of such separated condensates can act as a 
readout scheme for their relative phases [21] . We propose 
such a readout when one of the condensates is driven 
periodically with a frequency uj during the holdout for a 
single period T — 2tt/uj. Our specific prediction is that 
the relative phase measured for such a drive with uj = uj^ 
is going to match the phase without any drive indicating 
dynamic freezing. For all such experiments one needs 
to estimate a optimal temperature Tq at which they can 
be carried out. The typical value of U deep inside the 
Mott phase is ~ 2KHz = 200nK leading to a melting 
temperature of T™ ~ 0.2U = 40nK for d = 3. The SF 
phase near the Mott tip has a coherence temperature of 
Tc ~ zJc — 35nK [52]. Thus a temperature of a few 
nano-Kelvins (Tg ^ ^miPc): which is currently within 
the experimental reach, would be ideal for testing our 
prediction. 

In conclusion, we have demonstrated that periodic dy- 
namics of the ultracold bosons described by the Bose- 
Hubbard model leads to dynamic freezing of the Boson 
wavefunction at specific drive frequencies which are de- 
termined by the condition 4>d{T) = ATrm. The freezing 
phenomenon is qualitatively robust against the presence 
of the trap and quantum fiuctuations; it manifests itself 
at discrete drive frequencies < Jc via presence of dips 
in the defect density and can be detected by suitable in- 
terference experiments. 
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